Avalanches at rough surfaces 
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We describe the surface properties of a simple lattice model of a sandpile that includes evolving 
structural disorder. We present a dynamical scaling hypothesis for generic sandpile automata, and 
additionally explore the kinetic roughening of the sandpile surface, indicating its relationship with 
■ the sandpile evolution. Finally, we comment on the surprisingly good agreement found between 

this model, and a previous continuum model of sandpile dynamics, from the viewpoint of critical 
phenomena. 
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Avalanches are the signatures of instabilities on an evolving surface: regions on a sandpile, for example, which 
rotrude excessively from the surface, get dislodged by such mechanisms. This most intuitive picture of avalanching 
^ U is the one we seek to model and study: although presented here in the context of sandpiles, a similar picture may 
be relevant to intermittent granular flows along an inclined plane || , or to sediment consolidation Q . As deposition 
occurs on a sandpile surface, clusters of grains grow unevenly at different positions and roughness builds up until 
q ■ further deposition renders some of the clusters unstable. These then start 'toppling', so that grains from an already 
unstable cluster flow down the sandpile, knocking off grains from other similar clusters which they encounter. The net 
1— — '■ effect of this is to 'wipe off' protrusions (where there is a surfeit of grains at a cluster) and to 'fill in' dips, where the 
' oncoming avalanche can disburse some of its grains. In short, the surface is smoothed by the passage of the avalanche 
, so that there is a rough precursor surface, and a smoothed post-avalanche surface. 

In earlier work we have explored this issue via an analytical model involving coupled continuum equations. 

Here we use a cellular-automaton model @] of an evolving sandpile to look in more depth at the mechanisms by 
which a large avalanche smooths the surface. Our sandpile model is a 'disordered' and nonabelian version of the basic 
Kadanoff cellular automaton ||; a further degree of freedom, that involves granular reorganisation within columns, is 
added to the basic model which includes only granular flow between columns. In this sense, each column is regarded as 
' a cluster of grains so that we represent intracluster as well as intercluster grain relaxations, in accord with a previous 
"p-. , understanding of sandpile dynamics H . 

Our disordered model sandpile p0|| is built from rectangular lattice grains, that have aspect ratio a < 1 arranged 
in columns i with 1 < i < L, where L is the system size. Each grain is labelled by its column index i and by an 
orientational index or 1, corresponding respectively to whether the grain rests on its larger or smaller edge. 
The dynamics of our model have been described at length elsewhere but we review the essentials here: 

C ! 

• Grains are deposited on the sandpile with fixed probabilities of landing in the or 1 position. 



The incoming grains, as well as all the grains in the same column, can then 'flip' to the other orientation stochas- 
tically (with probabilities which decrease with depth from the surface). This 'flip', or change of orientation, is 
our simple representation of collective dynamics in granular clusters since typically clusters reorganise owing to 



the slight orientational movements of the grains within them |J . 

• Column heights are then computed as follows: the height of column i at time t, h(i,t), can be expressed in 
terms of the instantaneous numbers of and 1 grains, no(i,t) and ni(i,t) respectively: 

h(i, t) = ni(i, t) + ano(i, t) (1) 

• Finally, grains fall to the next column down the sandpile (maintaining their orientation as they do so) if the 
height difference exceeds a specified threshold in the normal way || (the pile is local, limited and has a fall 
number of two). At this point, avalanching occurs. 
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In earlier papers, ]7[, it was shown that the presence of disorder led, for large enough system sizes, to a preferred 
size of large avalanches. In the absence of disorder, scale-invariant avalanche statistics were observed. In this work, 
we focus on the state of the sandpile surface in a bid to correlate its evolution with the onset and propagation of 
avalanches. 



I. DYNAMICAL SCALING FOR SANDPILE CELLULAR AUTOMATA 

It is customary in the study of generalised surfaces to examine the widths generated by kinetic roughening [jl2| , and 
then establish properties related to dynamical scaling. However, the kinetic roughening of sandpile cellular automata 
has never been investigated; to begin with, therefore we postulate a principle of dynamical scaling for sandpile cellular 
automata in terms of the surface width W of the sandpile automaton: 

W(t) ~ t , t < t crossover = L z (2) 
W(L) ~ L a ,L -> oo (3) 

As in the case of interfacial widths, these equations signify the following sequence of roughening regimes: 

1. To start with, roughening occurs at the CA sandpile surface in a time-dependent way; after an initial transient, 
the width scales asymptotically with time t as t 13 , where (3 is the temporal roughening exponent. This regime is 
appropriate for all times less than the crossover time t crossover = L z , where z = a/P is the dynamical exponent 
and L the system size. 

2. After the surface has saturated, i.e. its width no longer grows with time, the spatial roughening characteristics 
of the mature interface can be measured in terms of a, an exponent characterising the dependence of the width 
on L. 

We define the surface width W(t) for a sandpile automaton in terms of the mean-squared deviations from a 
suitably defined mean surface; in analogy with the conventional counterpart for interface growth G^], we define the 
instantaneous mean surface of a sandpile automaton as the surface about which the sum of column height fluctuations 
vanishes. Clearly, in an evolving surface, this must be a function of time; hence all quantities in the following analysis 
will be presumed to be instantaneous. 

The mean slope < s(t) > defines expected column heights, h av (i,t), according to 

h av (i,t) = i < s(t) > (4) 
where we have assumed that column 1 is at the bottom of the pile. Column height deviations are defined by 

dh(i, t) = h(i, t) - h av (i, t) = h(i, t) - i < s(t) > (5) 
The mean slope must therefore satisfy 

Ei[h(i,t) -i < s{t) >] = (6) 

since the instantaneous deviations about it vanish; thus 

<s(t)>=2Ei[h(i,t)]/L(L + l) (7) 

(We note that this slope is distinct from the quantity < s'(t) >= h(L,t)/L that is obtained from the average of all 
the local slopes s(i,t) = h(i,t) — h(i — 1,<), about which slope fluctuations would vanish on average). 
The instantaneous width of the surface of a sandpile automaton, W(t), can be defined as: 



W(t) = y/Xi[dh(i,t) 2 ]/L (8) 

which can in turn be averaged over several realizations to give, < W >, the average surface width in the steady state. 
We also compute here the height-height correlation function, C(j,t), which is defined by 

C(j, t) =< dh{i, t)dh{i +j,t)>/< dh(i, tf > (9) 

where the mean values are evaluated over all pairs of surface sites separated by j lattice spacings: 

< dh(i,t)dh(i+j,t) >= T, i (dh(i,t)dh(i + j,t)) / '{L - j) (10) 

for < j < L. This function is symmetric and can be averaged over several realizations to give the average correlation 
function < C(j) >. 
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II. QUALITATIVE EFFECTS OF AVALANCHING ON SURFACES 



Before moving on to the quantitative descriptors of sandpile avalanching and surface roughening, we present some 
results using more qualitative indicators. Recent experiments || on sandpile avalanches have indicated that there 
are at least two broad categories; 'uphill' avalanches, which are typically large, and 'triangular' avalanches which 
are generally smaller in size. We have found evidence of this in a (2 + l)d disordered model of sandpile avalanches, 
which will be presented elsewhere fl3|| ; but in this work we discuss analogues in (1 + l)d, which are respectively 
'wedge-shaped' and 'flat' avalanches. The following data indicate that it is the larger wedge-shaped avalanches which 
alter surface slope and width, while the flatter, smaller avalanches alter neither very much. This is in accord with 
earlier work, where it was found that larger avalanches are the consequence of accumulated disorder, while the smaller 
ones can cause disordered regions to build up along the sandpile surface . 

Figure 0(a) shows a time series for the mass of a large (L — 256) evolving disordered sandpile automaton. The series 
has a typical quasiperiodicity [0]. The vertical line denotes the position of a particular 'large' event, while Figure 0(b) 




shows the avalanche size distribution for the sandpile. Note the peak, corresponding to the preferred large avalanches, 
which was analysed extensively in earlier work [[7|. Our data shows that the avalanche highlighted in Figure |l|(a) 
drained off approximately 5 per cent of the mass of the sandpile, placing it close to the 'second peak' of Figure 0(b). 
Figure 0(c) shows the outline of the full avalanche before and after this event with its initiation site marked by an 
arrow; we note that, as is often the case in one dimension, the avalanche is 'uphill'. The inset shows the relative 
motion of the surface during this event; we note that the signatures of smoothing by avalanches arc already evident as 
the precursor state in the inset is much rougher than the final state. Finally we show in Figure 0(d) the grain-by-grain 
picture of the aftermath pile superposed on the precursor pile, which is shown in shadow. An examination of the 
aftermath pile and the precursor pile shows that the propagation of the avalanche across the upper half of the pile 
has left only a very few disordered sites in its wake (i.e. the majority of the remaining sites are type) whereas the 
lower half (which was undisturbed by the avalanche) still contains many disordered, i.e. 1 type sites in the boundary 
layer. This leads us to suggest that the larger avalanches rid the boundary layer of its disorder-induced roughness, a 
fact that is borne out by our more quantitative investigations. 

In fact, our studies have revealed that the very largest avalanches, which are system-spanning, remove virtually all 
disordered sites from the surface layer; one is then left with a normal 'ordered' sandpile, where the avalanches have 
their usual scaling form for as long as it takes for a layer of disorder to build up. When the disordered layer reaches 
a critical size, another large event is unleashed; this is the underlying reason for the quasiperiodic form of the time 
series shown in Figure 0(a). 

Before moving on to more quantitative features, we show for comparison the sequence of Figure 0, for 

1. an ordered pile - Figure 0(a-d) 

2. a small disordered pile - Figure 0(a-d) 
We note the following features: 

• The small disordered pile has a mass time series (Figure 0(a)) that is midway between the scale-invariance of 
the ordered pile (Figure 0(a)) and the quasiperiodicity of the large disordered pile (Figure 0(a)). 

• The avalanche size distribution of the small disordered pile (Figure ||(b)) is likewise intermediate between that 
of the ordered pile (which shows the scale invariance observed by Kadanoff et al. (qj) and the two-peaked 
distribution characteristic of the disordered pile . 

• In both small and large disordered piles, we see evidence of large 'uphill' avalanches which shave off a thick 
boundary layer containing large numbers of disordered sites, and leave behind a largely ordered pile (see Figure 
0(c-d) and Figure 0(c-d)). By contrast the ordered pile loses typically two commensurate layers even in the 
largest avalanche, with a correspondingly unexciting aftermath state left behind in its wake (Figure 0(c-d)). 

We conclude from this that there is, even at a qualitative level, a post-avalanche smoothing of the sandpile surface, 
beyond a crossover length, as found in earlier work on continuum models H ; importantly, our discrete model reveals 
that this is achieved by the removal of (orientational) disorder, the implications of which we will discuss in our 
concluding section. The existence of the crossover length, in terms of the mass time series, has also been observed in 




experiment JlJ. 
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III. QUANTITATIVE EFFECTS OF AVALANCHING ON SURFACES 



A. Intrinsic properties of sandpile surfaces 

Inspired by the picture of smoothing avalanches, we have investigated many of the material properties of the sandpile 
in the special pre- and post-avalanche configurations. From these we have drawn the following conclusions: 

• The mean slope of the disordered sandpile peaks (see Table ||) before a large avalanche and drops immediately 
after; this statement is true for events of any size and thus remains trivially true for the ordered sandpile. 

• The packing fraction <f> of the disordered sandpile increases after a large event, i.e. effective consolidation occurs 
during avalanching (see Table |) . This consolidation via avalanching mirrors that which occurs when a sandpile 
is shaken with low-intensity vibrations fll5| , [T6| . 

• However, a far deeper statement can be made about the comparison of the surface width for pre- and post- 
large event sandpiles; Table | shows that the surface width goes down considerably during an event, once again 
suggesting that a rough precursor pile is smoothed by the propagation of a large avalanche. 

We have also investigated the dependence of various material properties of a disordered sandpile on the aspect ratio 
of the grains Table |l| shows our results, and Figure H illustrates the variation of the avalanche size distribution. 

There is a transition as aspect ratios of 0.7 are approached from above or below; we have shown above that piles 
with these 'critical' aspect ratios manifest strong disorder J?J in the sense of: 

• a 'second peak' in the avalanche size distribution denoting a preferred size of large avalanches 

• large surface widths denoting an increased surface roughness 

• a strong correlation between interfacial roughness and avalanche flow since the mean surface width varies dra- 
matically in the pre- and post- large event piles. 

Clearly, sandpiles containing grains with aspect ratios close to unity act essentially as totally ordered piles ||] ; there 
is however a significant symmetry in the shape of the avalanche size distribution curves above and below the transition 
region (see Figure ^(a) and (d)). These size distributions are reminiscent of those obtained in earlier work Q for the 
case of 'uniform disorder' (which referred to piles that have disorder throughout their volume rather than, as is the 
present case, disorder concentrated in a boundary layer). These observations lead us to speculate that there exist at 
least three types of avalanche spectra: 

1. the scale-invariant statistics characteristic of ordered sandpiles 

2. the strongly disordered statistics characterized by a second peak in the distribution; which we have obtained for 
specific values of the aspect ratio in the case where the disorder is concentrated in a boundary layer 

3. the more weakly disordered region (characterized by a flatter size distribution of avalanche sizes which is, 
nevertheless, not scale-invariant) obtained in the intermediate regimes of aspect ratio (as well as in the case of 
uniform disorder). 

It is clear that the presence of inherent inhomogeneities in grain shape ( which we describe quantitatively by aspect 
ratio) or bulk structure (which we describe by the classifications of 'uniform' or 'boundary' disorder) in a sandpile 
induces the presence of strong disorder in avalanche statistics. 

Additionally we present, in Figure ||, the mass-mass correlation function of a particular disordered sandpile; the 
curve has a peak, which indicates the average time between avalanches. Since the avalanche size distribution for this 
sandpile includes a preponderance of large events, we conclude that the peak in the correlation function corresponds 
approximately to the time between large avalanches. We also expect this timescale to manifest itself in the power 
spectrum of the avalanches; and we expect it to vary strongly with the level and nature of disorder in the sandpile. 
This work is in progress, as are efforts to relate the timescale found above to a characteristic spatial signature for 
large events. 

We present in Figure ^ the normalised equal-time height-height correlation function < dh(r + ro)dh(ro) > / < 
dh(ro) 2 > for a disordered sandpile. This shows that the height deviations (from the instantaneous expected column 
heights), in a disordered sandpile with L = 256, are positively correlated over about 80 columns, but also have a 
range where they are negatively correlated. In the inset we plot the related function 1— < dh(r + r )dh(r ) > / < 
dh(ro) 2 >=< (dh(r + ro) — dh(ro)) 2 > /2 < dh(ro) 2 >. For separations r much less than the correlation length of the 
system, we should have Jl8[ : 
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< (dh(r + r ) - dh(r )) 2 >~ |r| 2 ' 



(11) 



and, therefore, we would expect the function in the inset of Figure |6J to manifest a similar r-dependence. A linear fit 
to points with r < 30 (shown by the line in the inset of Figure ^ indicates a power-law dependence of the form r 67 
for r << L, implying that a ~ 0.34. As we will see below, this corresponds to the spatial roughening exponent of an 
ordered sandpile. An explanation of this behaviour is included in the next section. 



B. Spatial and temporal roughening of sandpile surfaces 

The hypothesis of dynamical scaling for sandpiles assumes that the roughening process occurs in two stages. First, 
the surface roughening is time-dependent, Eq. (Q); then once the roughness becomes temporally constant, the surface 
is said to saturate, and all further deposition results in surface fluctuations governed by Eq. (g). 

However, there is a subtlety concerning the first (i.e. time-dependent) stage; "early" model sandpiles are wedge- 
shaped and the transition to saturation is accompanied by a gradual build-up to a pile that has a single, sloping 
surface with a suitable angle of repose. 

We have taken this process into account to measure the dynamic exponent (3, Eq. (Erj; in this case surface widths 
are evaluated from the sloping portion of the pile. For the roughening exponent a, Eq. (]3|), we have measured surface 
widths from mature piles that have only a sloping surface. 

Our results are: 

• For disordered sandpiles (L = 2048) we find = 0.42 ± 0.05; for ordered sandpiles (L = 2048) (3 = 0.17 ± 0.05. 

• For disordered sandpiles above a crossover size of L c — 90 we find a = 0.723 ± 0.04; while for ordered piles we 
find a = 0.356 ± 0.05. 

• Based on the above values we find the dynamical exponent z, has values of 1.72 ± 0.29 and 2.09 ± 0.84 for the 
disordered and ordered sandpiles. 

The variation of the surface width, W, as a function of L, is shown in a log-log plot in Figure @. This fi gure shows 
clearly the crossover in a as a function of system size, for disordered sandpiles; the scaling behaviour of ordered 
sandpiles is shown for comparison. Disordered sandpiles with sizes below L c have a = 0.37 ± 0.05; this is in accord 
with earlier work, jl| where the second peak in the avalanche spectrum appeared only for disordered piles above 
crossover. The existence of this crossover length has been variously interpreted as a length related to reorganisation 
in the boundary layer of a sandpile || or to variations in the angle of repose in a (disordered) sandpile |hJ . Disorder 
appears to be crucial for the existence of such experimentally observed crossovers, since for example ordered models 
p0| , |2l| show no crossover in their measurements of a. The crossover effect in a disordered sandpile is also indicated 
by the height- height correlation function (Figure ^]). For separations r << L c ~ 90 in disordered sandpiles (with 
length L >> L c ), the exponent a obtained from the small r behaviour of the correlation function is that of the ordered 
sandpile. This suggests that even in disordered sandpiles grains which are within a crossover length, L Cl of each other 
tend to order i.e. the examination of the height-height correlation function for separations r << L c (Figure^) or the 
direct measurement of a for system sizes L << L c , as reported above, yields the exponent of the ordered sandpile, 
a -0.35. 

The above values indicate that while there is not a change of universality class as one goes from an ordered to a 
disordered sandpile (z stays the same, within the error bars), the disordered pile is clearly rougher with respect to 
both temporal and spatial fluctuations (a and f3 higher). 

It is important to note that our measurements of surface exponents are taken over many realisations of the surfaces 
concerned. Thus, even though, as demonstrated in earlier sections, the surface of a disordered sandpile is temporarily 
smoothed by the propagation of a large avalanche, it begins to roughen again as a result of deposition; the values of 
a and (3 that we measure are averages over millions of such cycles and hence reflect the roughening of the interface, 
in an average sense. By contrast, no abnormally large events occur for the ordered sandpiles and this is reflected by 
the lower values of fluctuations and exponents. 

The most striking aspect of these exponents is that they indicate that our present cellular-automaton model is a 
discrete version of earlier continuum equations J^], which were formulated independently, to model the pouring of 
grains onto a sloping surface. The exponents for our disordered pile are within error bars, exactly those that were 
measured for the height fluctuations of the surface in case 2, in ref. ||, while those for the ordered pile are exactly 
those that were measured for the fluctuations of the avalanches generated by the mobile grains, in the same case. 
This is in accord with the notion that the avalanches which flow on an ordered pile generate only mobile grains on 
the otherwise ordered surface, while as we have demonstrated above, avalanches that flow on a disordered pile, also 
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change the configuration of the surface by altering the distribution of height fluctuations (measured by the surface 
widths). We are exploring these analogies further, but note that this agreement is already a strong validation of both 
models. 



IV. DISCUSSION AND CONCLUSIONS 

We have presented a thorough investigation of the effects of avalanching on a sandpilc surface, focusing on the 
interrelationship between the nature of the avalanches and the surfaces they leave behind. We have also postulated 
a principle of dynamical scaling for sandpile surfaces, and measured the roughening exponents for a sample disor- 
dered sandpile. Finally, we have related the characteristics of avalanching in our model system to those obtained 
experimentally. 

Our current investigations concern several questions left unanswered above. These include the dependence of the 
crossover length L c on the disorder in the pile; as well as a fuller investigation of the effect of the nature of disorder 
(i.e. whether boundary or uniform). We would expect our correlation functions to depend strongly on the nature 
and magnitude of the disorder and we are undertaking a full quantitative study. Lastly, we hope that an extension of 
the present analysis to higher dimensions will yield more extensive comparisons with experiments than is presently 
available. 
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TABLE I. Instantaneous properties of disordered model sandpiles 



L 


State of pile 


Packing fraction (j> 


Slope 


Width 


256 


Before 


0.997 


1.15 


4.45 


256 


After 


1.000 


1.07 


2.06 


64 


Before 


0.991 


1.17 


1.41 


64 


After 


0.998 


1.04 


1.15 



TABLE II. Properties of model sandpiles 



Aspect ratio 


Packing fraction <j> 


Slope 


Width 


0.6 


0.997 


1.44 


2.33 


0.65 


0.997 


1.42 


2.34 


0.7 


0.997 


1.12 


3.76 


0.75 


0.997 


1.19 


3.74 


0.8 


0.998 


1.17 


2.35 


0.9 


0.999 


1.25 


2.29 


0.95 


1.000 


1.32 


2.37 


1.0 


1.000 


1.40 


2.41 
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FIG. 1. (a) A time series of the mass for a model sandpile (L = 256) that has been built to include a surface layer containing 
structural disorder (Mass is measured in lattice grain units). The vertical line indicates the position in this series of the 
large avalanche illustrated in Figure |l] c, d. (b) A log-log plot of the event size distribution for a model sandpile [L = 256) 
that includes a surface layer containing structural disorder, (c) An illustration of a large wedge shaped avalanche in a model 
sandpile (L = 256) that has been built to include a surface layer containing structural disorder. A lighter aftermath pile has 
been superposed onto the dark precursor pile and an arrow shows the point at which the event was initiated. The inset shows 
the relative positions of the two surfaces and their relationship to a pile that has a smooth slope (Height and position are 
measured in lattice units), (d) A detailed picture of the internal structure of a model sandpile in the aftermath of a large 
avalanche event. The individual grains of the aftermath pile (for columns 1 — 128 of a sandpile with L = 256) are superposed 
on the gray outline of the precursor pile. 

FIG. 2. (a) A time series of the mass for a model sandpile (L = 256) that has been built to exclude any structural disorder 
(Mass is measured in lattice grain units). The vertical line indicates the position in the series of the avalanche illustrated 
in Figure ^ c, d. (b) A log-log plot of the event size distribution for a model sandpile (L — 256) that excludes structural 
disorder, (c) An illustration of a large 'flat' avalanche in a model sandpile (L = 256) that excludes structural disorder. A 
lighter aftermath pile has been superposed onto the dark precursor pile and an arrow shows the point at which the event was 
initiated. The inset shows the relative positions of the two surfaces and their relationship to a pile that has a single smooth 
slope (Height and position are measured in lattice units), (d) A detailed picture of the internal structure of a model sandpile 
in the aftermath of a large avalanche event. The individual grains for the aftermath pile (for columns 1 — 128 of an ordered 
pile with L — 256) are superposed on the grey outline of the precursor pile. 
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FIG. 3. (a) A time series of the mass for a model sandpile (L = 64) that has been built to include a surface layer containing 
structural disorder (Mass is measured in lattice grain units). The vertical line indicates the position in this series of the large 
avalanche illustrated in Figure ^ c, d. (b) A log-log plot of the event size distribution for a model sandpile (L = 64) that includes 
a surface layer containing structural disorder, (c) An illustration of a large wedge shaped avalanche in a model sandpile (L = 64) 
that has been built to include a surface layer containing structural disorder. A lighter aftermath pile has been superposed onto 
the dark precursor pile and an arrow shows the point at which the event was initiated. The inset shows the relative positions 
of the two surfaces and their relationship to a pile that has a smooth slope (Height and position are measured in lattice units) . 
(d) A detailed picture of the internal structure of a model sandpile in the aftermath of a large avalanche event. The individual 
grains of the aftermath pile are superposed on the gray outline of the precursor pile. 

FIG. 4. Exit mass event size 

distributions for a disordered sandpile model with L = 256, and a — a) 0.6, b) 0.7, c) 0.75, d) 0.85, e) 1.0. The curves have 
each been shifted by 0.5 to make them distinct. 



FIG. 5. The mass-mass correlation function for a disordered model sandpile with L — 256. 

FIG. 6. The normalized correlation function of column height deviations for a disordered model sandpile with L — 256. The 
inset illustrates an initial decay with r, going as r -0 ' 67 (r is measured in lattice units). 



FIG. 7. A log-log plot of the surface widths, W, against the system size, L, for model sandpiles; (x) ordered piles, (o) 
disordered piles. Widths are measured in lattice units and all points have error bars that are ~ 0.02. 
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